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In this work, the spin relaxation accompanying the spin diffusion in symmetric Si/SiGe quantum 
wells without the D'yakonov-Perel' spin-relaxation mechanism is calculated from a fully microscopic 
approach. The spin relaxation is caused by the inhomogeneous broadening from the momentum- 
dependent spin precessions in spatial domain under a magnetic field in the Voigt configuration. 
In fact, this inhomogeneous broadening together with the scattering lead to an irreversible spin 
relaxation along the spin diffusion. The effects of scattering, magnetic field and electron density on 
spin diffusion are investigated. Unlike the case of spin diffusion in the system with the D'yakonov- 
Perel' spin-orbit coupling such as GaAs quantum wells where the scattering can either enhance 
or reduce spin diffusion depending on whether the system is in strong or weak scattering limit, 
the scattering in the present system has no counter-effect on the inhomogeneous broadening and 
suppresses the spin diffusion monotonically. The increase of magnetic field reduces the spin diffusion, 
while the increase of electron density enhances the spin diffusion when the electrons are degenerate 
but has marginal effect when the electrons are nondegenerate. 

PACS numbers: 72.25.Rb, 72.25.Dc, 71.10.-w 



I. INTRODUCTION 

The study of semiconductor spintronics has attracted 
a great deal of attention for its potential application 
to spin-based devices.— Among different prerequi- 
sites for realizing these devices, such as efficient spin 
injection^i^ and suitable spin lifetime^si long spin dif- 
fusion/transport lengtb^'^'i^ is required sometimes, es- 
pecially for the design of spin transistor and spin valve. 
Therefore, it is important to investigate the spin diffu- 
sion/transport in semiconductors. 

In the study of spin diffusion/transport, the two- 
component drift-diffusion model is widely used in the 
literaturci^ i^^'^^i-'^^'^'^i^^'^^ In this model the spin-diffusion 
length Ls is connected to spin-relaxation time Tg through 
spin-diffusion coefficient 0^: Lg — ^/]J~Ts ^ilLlLll 
with Ds usually assumed to be equal to the charge- 
diffusion coefficient Z?j,^ii^iiiii^ii^ This equation im- 
plies infinitely long spin-diffusion length when the spin- 
relaxation time Ts goes to infinity. In bulk Si, there is no 
D'yakonov-Perel' (DP) spin-orbit couplingil. due to the 
bulk inversion-symmetry. Thus the DP spin-relaxation 
mechanismii is absent and the spin-relaxation time in 
bulk Si is infinite when the other spin-relaxation mech- 
anisms (such as the Elliott- Yafet mechanisroi^) are ig- 
nored. Therefore, an extremely long or even infinite spin- 
diffusion length is expected in bulk Si. However, this is 
not the case in the presence of a magnetic field. Very 
recently, Appelbaum et al. studied the spin transport 
in bulk Si with a magnetic field perpendicular to both 
directions of spin transport and spin polarizatioUfi^ It 
has been shown that the spin-diffusion length is very 



small when the magnetic field becomes slightly strong 
(typically the spin-diffusion length is about several mi- 
crons when the magnetic field is in the order of 0.1 T)J^ 
To account for the experimental spin relaxation and de- 
phasing (R&D) along spin transport as well as the small 
spin-diffusion length, the drift-diffusion model was uti- 
lized and the interference among different spin-precession 
angles when reaching the same distance due to dis- 
tinct transit times was suggested to be important 
In fact, this spin R&D along spin diffusion due to the 
magnetic field in the absence of the DP spin-relaxation 
mechanism was predicted from a fully microscopic ap- 
proach, i.e., the kinetic spin Bloch equation (KSBE) 
approach^>S2^^^^ back in 2002^ In this ap- 
proach, the momentum-dependent spin precessions give 
rise to the inhomogeneous broadeningi22i2i In the pres- 
ence of the inhomogeneous broadening, any scattering 
(including the spin-conserving scattering) leads to an ir- 
reversible spin R&Di^i2iiS^ In spin diffusion/transport, 
it has been shown that the inhomogeneous broadening 
is determined by the spin-precession frequency lO]^ = 
■p^(r2k + gi^B^) when the spin diffusion/transport is 
along the x-axisjSS Here fJk is the DP term.— It was 
shown in Ref. [13 that even when ilk = 0, the kx de- 
pendence in u;k still causes spin R&D. This is exactly 
the case in the experimental work on Si. The present 
work is to investigate the spin diffusion in Si/SiGe quan- 
tum wells (QWs) by means of the KSBE approach, in 
order to gain a deeper insight into the spin relaxation 
along spin diffusion in the absence of the DP spin- 
relaxation mechanism. The QWs are symmetric with 
even number of monoatomic Si layers and ideal het- 
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erointerfaces to exclude the Rashba spin-orbit coupling. 
In addition, the study of spin relaxation in asymmetric 
Si/SiGe QWs has been carried out theoretically^^ and 
experimentally,™i22, showing that the Rashba spin-orbit 
coupling"^^ is very small (typically about three orders of 
magnitude smaller than that in QW structures based on 
III-V semiconductors^S) and the spin-relaxation time is 
quite long (in the order of lO"'^ - lO'^ s),^^^^^ Thus 
even for asymmetric Si/SiGe QWs, the present study still 
makes senses as long as the DP term fJk is weak enough 
compared to the magnetic-field term gus^- 



II. MODEL AND KSBES 

We start our investigation from an n-type symmetric 
Si/SiGe QW with its growth direction along the 2||[001] 
direction. The lowest conduction band in bulk Si is lo- 
cated near the X points of the Brillouin zone. Due to 
the quantum confinement along the z-direction, the two 
degenerate Xz valleys lie lower than X^ and Xy valleys. 
The well width a is set as small as 5 nm and the tempera- 
ture is lower than 80 K, thus with the moderate electron 
concentration, only the lowest subband of X^ valley is 
occupied. The spin polarization is injected constantly 
from one side of the sample {x — Q plane) with polar- 
ization P and diffuses along the x-axis, while the spatial 
distribution of spin in the y-direction is uniform. The 
magnetic field B is applied in the x-y plane. It can be 
along arbitrary direction without inducing any essential 
difference in the absence of the DP spin-orbit coupling. 
However, for the sake of convenience, the magnetic field 
is set to be along the y-axis. 

As the two Xz valleys are degenerate, only one of them 
needs to be considered. However, both the intra- and 
inter- valley scatterings have to be taken into account. 
The KSBEs for one valley read2^i^i25i26 
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Here pis^{x,t) represent the density matrices of elec- 
trons with two-dimensional momentum k (referring to 
the bottom of Xz valley under investigation) at posi- 
tion X and time t. Their diagonal terms Pk,(T<j=/k,cr 
{a— ±1/2) represent the electron-distribution functions 
and the off-diagonal ones /3k,i/2,-i/2=Pk -1/2 1/2 
scribe the inter-spin-band correlations for the spin co- 
herence. Ek(a;,i) is the Hartree-Fock term from the 
Coulomb interactioni^i^ii^ — e is the electron charge, 
mt = 0.196too is the transverse effective mass in the x-y 
plane, and g « 2 is the effective g factor for electrons 
in X valleys of Sii^^ '^{x,t) is the electric potential de- 
termined by the Poisson equation ^ g^f'*^ = e[n{x,t) — 



A''o]/(47reoKoa) with n{x,t) = ^^na{x,t) standing for 
the electron density at position x and time t. Nq is the 
background positive charge density, and n{x,0) = Nq 
denoting the initial uniform spatial distribution of elec- 
tron density, kq = 11.9 is the relative static dielectric 

, intra ^ , ^, inter 

and gt' originate from 
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the intra- and inter-valley scatterings, respectively. They 
are composed of the electron-phonon, electron-impurity 
and electron-electron Coulomb scatterings (see also the 
Appendix) . 

Before numerically solving the KSBEs, a much sim- 
plified situation with the elastic electron-impurity scat- 
tering only is investigated analytically. Based on this 
simplified investigation, some properties of spin diffu- 
sion in Si/SiGe QWs can be speculated. Assuming 
k = fc(cos 9, sin 6) , in the steady state the Fourier compo- 
nents of the density matrix with respect to angle 9 obey 
the following equation: 
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Here 



J^^ d9{l — cosl9)U^ is the rth-order momentum- 
relaxation rate, with |q| = y/2k'^{l — cos9). is the 
electron-impurity scattering potential and Ni represents 
the impurity density. It is noted that ~ and 

^ = 0. From the above equation one can obtain a closed 

group of first-order differential equations for p^^ and p^ 
by neglecting higher orders of p'^, with \l\ > 1. From these 
equations the following second-order differential equation 
about p1 is obtained: 

0/ N of mtgp-BBV 



\hk 



(3) 



Defining the "spin vector" as S°(a:) = Tr[p^{x)(T\ and 
using the boundary conditions (i) S5J.(0) = (0,0, S*?)^ 
^{ 

k\ 



and (ii) S"^,(-t-cx) 



0, S°(x) is found to be 



/ . 
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with ujp = gpsB/h being the spin-precession frequency 
in time domain under magnetic field B = By. The to- 
tal spin signal in the z-direction (including two valleys) 
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reads S^x) = 2 Tr[pk(x)a,] = J+°° fhS^^ix). To 
get 52(2;), S*^^ in boundary condition (i) is assumed to 

be = fk,l/2 - /fe -1/2 with fk,cr = ^mn^k'^ 

Here /i^ is determined by -1/2) = Nq 

and /+°°^ f fc(/fe4/2 - fK-1/2) - ^A^o with P the spin 
polarization at a; = plane. 

It is noted from Eq. ([4]) that during the spin diffu- 
sion, S^(a;) on one hand precesses around the direction 
of magnetic field, and on the other hand decays in magni- 
tude. However, when becomes infinity, i.e., in the limit 
of zero electron-impurity scattering, Eq. ^ becomes 

Slix) = Si (sin ^^.,0, cos ^^x)''. This solu- 

tion clearly indicates the momentum dependence of spin 
precession in spatial domain under magnetic field, which 
leads to the inhomogeneous broadening. ^'^ This inhomo- 
geneous broadening alone leads to a reversible decay of 
the total spin signal Sz{x) along the x-direction. How- 
ever, in the presence of scattering, this decay becomes 
irreversible, as shown in the exponential damping term 
for each Sl^{x) in Eq. Q. Therefore, even without the 
DP term, there is spin relaxation accompanying the spin 
diffusion. Due to the factor S^^ included in Eq. ([4]), the 
main contribution to Sz {x) comes from fc-states near the 
Fermi surface at x = plane, especially when the tem- 
perature is low. Thus 



indicating that the spin diffusion can be suppressed by 
the scattering strength and the magnetic field, based on 
the increase of the exponential damping rate with I/t^^. 
and ujp. It has been shown earlier that, in systems with 
the DP spin-orbit coupling, the scattering on one hand 
provides the spin-R&D channel in the presence of the 
inhomogeneous broadening from the DP term but, on 
the other hand, can have a counter-effect on the inho- 
mogeneous broadening 3ii22iS£ In the strong scattering 
limit, the counter-effect dominates and thus increasing 
scattering strength results in a suppression of spin R&D 
in time domain^iiS^ or an enhancement of spin diffu- 
sion/transport in spatial domain. However, for the cur- 
rent case without the DP spin-orbit coupling, the scatter- 
ing shows no counter-effect on the inhomogeneous broad- 
ening but just suppresses the spin diffusion. Except for 
the above analysis based on the simplified model, to gain 
a complete picture of the problem of spin diffusion, one 
must solve the KSBEs taking into account the electron- 
phonon and electron-electron Coulomb scatterings, both 
of which play an important role on spin R&D. 



III. NUMERICAL RESULTS 

To numerically solve the KSBEs, the double-side 
boundary conditions are used.— These conditions assume 
a steady spin polarization P for electrons with fe^; > at 
boundary x — and a vanishing spin polarization at fi- 
nite sample length x — L for electrons with kx < 0^ For 
the Poisson equation, the boundary conditions are set to 
be ^(0,t) = ^{L,t) = 0. The numerical scheme for solv- 
ing the KSBEs is given in Ref. 26. Once the KSBEs are 
numerically solved, the steady-state distribution of elec- 
trons Na-{x) = no- (a;, +00) = 2J2-kf'k.a-{x,+oo) (factor 2 
comes from the two degenerate Xz valleys) is calculated 
and thus the spin signal Sz{x) — Ni/2{x) — N_i/2(x) 
can be derived to analyze the spin-diffusion properties. 
In the calculation, the electron density is set to be 
4x10^^ cm~^ except otherwise specified, and the impu- 
rity density Ni is assumed to be O.ITVq when the impu- 
rities are present. Furthermore, the initial spin polariza- 
tion P at X = plane is set to be 5 %. The effects of 
scattering, magnetic field and electron density on spin 
diffusion are investigated, with the main results given in 
Figs.IlEl 

We first investigate the effects of scattering and mag- 
netic field on spin diffusion. The steady-state spatial dis- 
tributions of spin signal calculated with different scat- 
tering are shown in Fig. [TJ In order to show the prop- 
erty of spin diffusion clearly, we also plot the absolute 
value of Sz vs. a; on a log-scale in the same figure. 
All these curves indicate obvious spin relaxation along 
spin diffusion without the DP spin-relaxation mecha- 
nism. By comparing the curves labeled as "EE" and 
"EI" , one finds that the electron-electron Coulomb scat- 
tering can suppress spin diffusion effectively. In fact, 
the Coulomb scattering plays an important role in both 
spin R&D2ii2^i2i and spin diffusion/transport.— s^i^Si In 
a system with the DP spin-orbit coupling, the Coulomb 
scattering not only contributes to the total momentum 
relaxation time t^,'^^ but also has a counter-effect to 
the inhomogeneous broadeningi ^-'^i^^i^'^'^^i^^ Therefore, in 
the strong scattering limit, adding Coulomb scattering 
may suppress the spin relaxation and enhance the spin 
diffusion/transport For the current situation with- 
out the DP spin-orbit coupling, the Coulomb scatter- 
ing affects the spin diffusion only through and thus 
only suppresses the spin diffusion/ transport. Similarly, 
the electron-phonon scatterings, including both the intra- 
and inter- valley scatterings, also contribute to the mo- 
mentum relaxation and suppress spin diffusion effectively, 
as shown by the "EP" curve. For the case of stronger 
scattering strength (i.e., with all the different scatter- 
ings included), the spin-diffusion length becomes much 
smaller. It is also noted that with the increase of scatter- 
ing strength, the spatial spin-precession period decreases. 
This is because the spin-precession frequency increases 
with as shown in Eq. ([5]). The magnetic- field de- 

pendence of spin diffusion is investigated as well, with the 
results corresponding to different magnetic fields shown 
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in Fig.[2l It is revealed that both the spin-difFusion length 
and the spin-precession period decrease with an increase 
of magnetic field strength B. This is because both the 
damping rate and the spin-precession frequency increase 
with ujp and thus B, as shown in Eq. ([5]). These studies 
indicate a suppression of spin diffusion due to the scat- 
tering and magnetic field, just as obtained earlier with 
the simplified model. 
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FIG. 1: (Color online) The steady-state spatial distributions 
of spin signal Sz with different scatterings included. The 
curve labeled "EE", "EP" or "EI" stands for the calcula- 
tions with the electron-electron, electron-phonon or electron- 
impurity scattering, respectively, while the curve labeled 
"El-f EE-f EP" stands for the calculation with all the scatter- 
ings. In order to get a clear view of the decay and precession 
of Sz, we also plot the corresponding absolute value of Sz 
against a:: on a log-scale (Note the scale is on the right hand 
of the frame) . The dashed curves correspond to the part with 
Sz < 0. 

The density dependence of spin diffusion is also inves- 
tigated. The steady-state spatial distributions of spin 
signal with three different electron densities under tem- 
perature T = 40 and 10 K are plotted in Fig. [3)Ja) and 
(b) , respectively. The spin signal is rescaled by the corre- 
sponding electron density to be Sz/Nq for comparison. It 
is noted that when T = 40 K, the spin diffusion is insensi- 
tive to the electron density. That is because the electrons 
are nondegenerate in the studied electron-density regime 
when T = 40 K, due to the large transverse effective elec- 
tron mass in X valleys of Si. However, when T decreases 
to 10 K, the effect of electron density on the spin dif- 
fusion can be seen, as shown in Fig. [3Jb). In the large 
density regime where the electrons become degenerate, 
the spin diffusion length increases with the electron den- 
sity (compare the situations with Nq = 4.0 x 10^^ and 
l.Ox 10^^ cm^^). This is mainly due to the decrease of the 
damping rate with kf, as shown in Eq. In addition, 
in the low density regime it is shown that the density 
again has a marginal effect on spin diffusion (compare 
the situations with Nq = 1.0 x 10" and 0.5x10" cm'^), 
as the electrons remain nondegenerate there. It is noted 
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FIG. 2: (Color online) Sz vs. x in the steady state for different 
magnetic field strengths. Solid curve: B = 2T; Dashed curve: 
B = 1 T; Dotted curve: B = 0.5 T. T = 80 K and Ni = 
O.liVo. 



that the density dependence of spin diffusion in Si is very 
different from the density dependence of the spin relax- 
ation in GaAs QWs where non-monotonic density depen- 
dence was predicted^^ and realized experimentally very 
recently."^^ This is due to the fact that in GaAs QWs 
the inhomogcncous broadening comes from the DP term 
which is, however, absent in the symmetric Si/SiGe QWs. 
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FIG. 3: (Color online) Sz/Nq vs. x in the steady state with 
different electron densities at (a) T = 40 K and (b) T = 10 K. 
B = 1 T and Af, = 0. 
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IV. SUMMARY 

In summary, the present work investigates the spin dif- 
fusion in symmetric Si/SiGe (001) QWs at low tempera- 
ture. There is no DP spin-relaxation mechanism due to 
the absence of the DP spin-orbit coupling in this system. 
However, a magnetic field in the Voigt configuration is 
present. Our simulations were performed in a fully mi- 
croscopic way based on the KSBE approach, with all the 
relevant scatterings included. It was shown that, even 
without the DP spin-relaxation mechanism, the electron 
spins relax effectively along the spin diffusion. This spin 
relaxation is caused by the inhomogeneous broadening 
from the momentum-dependent spin precessions in spa- 
tial domain. The effects of scattering, magnetic field and 
electron density on spin diffusion were investigated. It 
was shown that, unlike the case of spin diffusion in the 
system with the DP spin-orbit couplingj^^ in Si/SiGe 
(001) QWs any scattering suppresses the spin diffusion 
without any counter-effect on the inhomogeneous broad- 
ening. The magnetic field reduces spin diffusion also. It 
was further revealed that the increase of electron density 
enhances the spin diffusion when the electrons are de- 
generate but has marginal effect when the electrons are 
nondegenerate. 
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APPENDIX A: THE SCATTERING TERMS OF 
THE KSBES 



The scattering terms are analogous to those shown 
in Ref. iSfi, except the following differences. The two 
valleys are degenerate here, thus the multivalley KS- 
BEs similar to those shown in Ref. [H can be sim- 
plified to obtain the KSBEs of a "single" valley [i.e., 
Eq. ([T])]. The intra- valley scatterings consist of those due 
to longitudinal-acoustic (LA) and transverse-acoustic 
(TA) phonons, while the inter-valley scatterings are the 
g-type scatterings involving LA, TA and longitudinal- 
optical (LO) phonon branches.'^'' However, due to the 
low temperature, the scattering due to the LO-phonon 
is neglected. M^ ^,^^ ^ = _g£_|/.„j^^(ig^)|2 jg 
the matrix element for the intra-valley scattering, and 

^a,mter,Q = 2da„.,„"e..Q l-^"*^'' ) 1^ ^^'^ iutcr-valley 
scattering, a — LA/TA stands for the LA/TA phonon 
mode, d = 2.33 g/cm"^ is the mass density of Sii^ 
Dla = 6.39 eV and Dta = 3.01 cY.^l r2a,mtra,Q = VaQ 
with phonon velocities vla — 9.01 x 10^ cm/s and 
VTA = 5.23 X IQS cm/s.37 Ala = 1-5 x 10^ eV/cm 
and At A = 0.3 x 10® eV/cm.^^ The phonon energies for 
the inter-valley scattering are approximately fixed to be 

fl^LA,inter,Q = 0.019 cV and hVLTA,inter,Q. = 0.01 gV ML. 

\I~f{iqz)\'^ — y^^yf^Jiyi (7 — intra/inter) is the form fac- 
tor with y = aqz/2 for the intra-valley scattering and 
y = a{qz — 2K^J/2 for the inter- valley scattering. Here 
= 0.85 X ^ with uq being the Si lattice constant is 
the z component of the coordinate of the bottom of 
valley. 
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